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Abstract 

The quantized vortex state is investigated in a Bose-Einstein condensate, 
confined in a multiply connected geometry formed by a Laguerre-Gaussian 
optical trap. Solving the Gross-Pitaevskii equation variationally, we show that 
the criterium for vortex stability is that the interatomic interaction strength 
must exceed a critical interaction strength. The time evolution of a freely 
expanding Laguerre-Gaussian condensate with a vortex is calculated, and 
used to derive the interference pattern of such a condensate overlapping with 
a parabolically trapped condensate. 
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I. INTRODUCTION 



Superfluidity, and its characteristic manifestation as a state with quantized circulation 
(a vortex), is intimately connected with the phenomenon of Bose-Einstein condensation, yet 
the precise relation between superfluid persistent currents and Bose-Einstein condensation 
(BEC) is only beginning to be elucidated PJ. Before 1995, the link between superfluidity 
and BEC was almost exclusively studied in the context of liquid helium ( 4 He and 3 He) 
0, where the study of the relation between superfluidity and BEC is complicated by the 
strong interatomic interactions between the atoms in the liquid. In 1995, Bose-Einstein 
condensation was realized in magnetically trapped clouds of alkali atoms |§. In these novel 
condensates, the bosonic atoms are weakly interacting (contrary to the case of liquid helium), 
and as such these systems have the potential to shed new light on superfluidity. Soon after 
the initial creation of alkali gas condensates, several experimental groups set out to create 
a vortex - a quantum of superfluid circulation - in this novel system Q. 

The initial attempts to create a vortex in a condensate by stirring the trapped condensate 
with a blue detuned laser beam || were unsuccessful, even though early theoretical work by 
Dalfovo et al. indicated that persistent superfluid currents can indeed manifest themselves 
in Bose-Einstein condensates as vortices 0. Subsequent analysis || showed that vortices 
are unstable in the simply-connected, not-stirred, spin-polarized condensates formed in the 
original experiments |§. 

This can be understood as follows. Along the vortex line, the order parameter of the 
Bose-Einstein condensate has to vanish. Phrased metaphorically, a "hole" has to be "drilled" 
in the condensate along the vortex line. In the magnetic trap, the modulus square of the 
order parameter of a condensate without a vortex is largest in the center of the trap. As a 
consequence, a vortex line through the center of the trap will perturb the order parameter 
more than a vortex line at the edge of the condensate (since, using our metaphor again, it 
will cost more energy to "drill" through the center of the condensate than through its edge). 
It will be energetically favorable for the vortex line to be at the edge of the condensate. 
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In the presence of dissipation, this will cause the vortex line to migrate to the edge of the 
condensate so that the vortex condensate will decay into a non-vortex state. This argument, 
sketched here with some roughness, has been worked out with precision by Fetter and co- 
workers in [^]. 

Several schemes have been proposed, both theoretically and experimentally, to stabilize 
vortices: rotating the trapping potential (analogous to rotating a bucket containing 4 He) 
H or stirring the condensate with an off- resonance laser ||, raising the temperature (to 
'pin' the vortex in the potential created by the non-condensate fraction at the center of 
the vortex) ||, phase imprinting methods flOfl , and various other techniques ||11|| . Recently, 
vortices were created experimentally, both with a 'rotating bucket' experiment El and with 



the use of a two-component condensate |13 |. In the latter experiment, one of the components 
of the spinor condensate 'pins' the vortex present in the other component and a Ramsey type 
interference between the two components is used to detect the vortex. The long lifetimes 
of the two-component condensates in ref. |13[ is due to similar singlet and triplet scattering 
lengths of rubidium, resulting in an anomalously low inelastic loss rate [T^ . 



Nevertheless, stable vortices have not yet been realized in a spin-polarized condensate in 
non-rotating traps, and new methods of stabilization and detection must be developed. A 
promising scheme, based on an analysis similar to that of Fetter f7j, is the use of multiply- 
connected condensates. A candidate trap to create a multiply connected condensates consists 
of a red-detuned laser beam in a Laguerre- Gauss mode |L5],|16|, which we discuss in section 
II. If the condensate order parameter has a toroidal geometry, a vortex line threaded through 
the cylindrical symmetry axis of the torus will not perturb the condensate order parameter 
strongly, and moreover create a metastability barrier for vortex decay. One of the goals of 
this paper is to verify this statement about vortex stability (in section III). In a condensate 
with a toroidal order parameter, the vortex line can no longer be detected as a line along 
which the density of Bose-Einstein condensed atoms vanishes. To detect vorticity, we propose 
(in section IV) a method based on interference, similar to that proposed in refs. f° r 
simply connected condensates. 



II. BOSE-EINSTEIN CONDENSATION IN LAGUERRE GAUSS TRAPS 



A. Trapping geometry 



Toroidal confinement for ultracold atoms can be obtained by an optical dipole trap flS 
which consists of a laser beam in a Laguerre-Gaussian mode. The Laguerre-Gaussian mode 
{n, m} is characterized by an intensity profile given by [EC 
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where r is the radial distance from the center of the beam, z is the position along the 
propagation direction of the beam, Wq (referred to as the 'waist parameter') is a parameter 
controlling the minimal width of the beam, k is the wave number of the laser, and L n j is 
the Laguerre polynomial of order {n, I}. Such laser beams have modes which show a node 
in the center, and trap the atoms in a cylindrical shell around the axis of propagation of the 
beam. The intensity profile of the laser beam in the Laguerre-Gaussian mode is illustrated 
in Figure 1. Laguerre-Gaussian (LG) laser beams have already been used to successfully 



trap atoms [16], and subsequent theoretical work has shown that toroidal traps formed by a 



red-detuned LG beam can be loaded from initial conditions similar to those of conventional 



magnetic traps [21] . Once the trap has been loaded, one possibility to create a vortex in the 
trapped gas would be by a phase imprinting method, which already successfully resulted in 
the creation of solitons in condensates WM- Alternatively, a rotating perturbing potential 



can be used to stir the condensate and set up persistent flow in the toroidal geometry 

In this paper we investigate the properties of the vortex Bose-Einstein condensate, op- 
tically trapped by a laser beam in a Laguerre-Gaussian propagation mode {n, 1} = {0, 1}. 
Along the z-axis (the direction of propagation of the laser beam) an additional magnetic trap 



15] results in a parabolic z-axis confinement with frequency f2, which prevents the atoms 
from escaping along the direction of propagation of the laser beam. This 'plugging' of the 
optical trap was achieved in [|16[ by using blue-detuned 'plugging beams'. The condensate in 
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the Laguerre- Gauss geometry |15| will be denoted by 'Laguerre-Gaussian condensate' (LG 



condensate), in contrast with the condensate in a parabolic confinement. 

The intensity profile of such a Laguerre- Gauss beam in the zy-plane (where z is the axis 
of propagation of the laser beam), given by ([I]), is shown in Figure 1. In the remainder of 
this paper, we use units so that h = m = Q = 1 (where m is the mass of the atoms). In 
these units, the laser beam parameters for Figure 1 are chosen as follows: k = 2, Wo = 5. 
Figure 2 illustrates a surface of constant intensity of the laser beam; a hollow cylindrical 
shell. This will also be the shape of the cloud of trapped atoms The extension of the 



cloud along the z-axis can be tuned by selecting the Q frequency of the magnetic trap along 
the z-axis. 

B. Mean- field approach 

Confined Bose-Einstein condensates are well described by a mean-field theory where the 
properties of the condensate are derived from a complex function This function (the order 
parameter) is interpreted as a macroscopic wave function and obeys the Gross-Pit aevskii 
equation f24| : 

_ |!_A* + Vconf (p, + U \^ = EV. (2) 
2m 

The term nonlinear in \I/ in (fj) arises from the interparticle interaction potential, which is 
treated as a contact potential with scattering length a scat so that Uq = 47rh 2 a scat /m where 
m is the mass of an atom. The optical confinement is generated by a red-detuned Laguerre- 
Gaussian laser beam in the {0, 1} mode, with waist parameter Wq and wave number k, and 
propagating along the z-axis |[20|| . In addition to this, a harmonic confinement (characterized 



by a frequency Q) is present which confines the atoms along in the ^-direction to a region 
z < kWl [0. The potential energy corresponding to this geometry is given, in cylindrical 
coordinates {p,9,z}, by: 

WP,M) = -A ^exp{-2p 2 /^ 2 } + ^z\ (3) 



In the present treatment, the Gross-Pitaevskii equation (Q) will be solved variationally J|5 . 



Consider for this purpose the following two trial functions, expressed in cylindrical coordi- 
nates {p, 9, z}\ 

VO, 9, z) = Af exp{-a (p - Po) 2 /2 - Co^ 2 /2}, (4) 
^ v (p,9,z)=Af 1 pexp{- t r 1 (p- Pl ) 2 /2 - ( lZ 2 /2}e w . (5) 

The variational parameters in the trial functions are 00,01, controlling the width of the 
functions in the radial direction; CqjCi> controlling the widths in the ^-direction; and Po,Pi, 
controlling the radial displacement of the maximum of the function away from the z-axis. 
Ao and A/"i are normalization constants, which are determined by demanding that the norm 
of the trial function equals the number of particles in the condensate. 

The circulation around a closed loop in the trapped, interacting Bose gas is defined 
by k = ^§ VS 1 , where 5* represents the phase of the order parameter solving the Gross- 
Pitaevskii equation (|2|). The single- valuedness of the order parameter ensures that the 
circulation in the trapped Bose gas is quantized: n = nh/m with n = 0, 1, 2, .... A vortex is 
present whenever n > 0. The trial function ip(p, 9, z) has a constant phase and hence there 
is no circulation, no vortex. On the other hand, the phase of ip v (p, 9, z) changes by 2ir along 
any closed loop encircling the z-axis. Thus ip v is a trial function for the order parameter 
of a condensate with one quantum of vorticity: its circulation equals k = h/m. The trial 
function i[) v must have a node along the z-axis (the factor p in expression (||)): if this node 
were not present, the variational kinetic energy would diverge along the z-axis. 

We have used ip and ip v as variational trial functions to solve the Gross-Pitaevskii equa- 
tion for a condensate without and with a vortex, respectively. The resulting variational 
energy is shown in Figure 3 as a function of Na scat /(iHo where N is the number of atoms 
and ciho — a/ h/ (rrdfj. The quantity Na scat /a HO is a dimensionless measure of the strength 
of the interaction: upon increasing Na scat / 'auo > the interatomic interaction becomes 
more repulsive. The energy of the LG condensate with a vortex is larger than the energy of 
the LG condensate without a vortex for all interaction strengths investigated. The difference 
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is smallest for the non- interacting gas and increases monotonically as Na sca t/ clho increases. 
This means that for all investigated interactions strengths the vortex LG condensate is ei- 
ther metastable or unstable. In inset (a) of Figure 3, the optimal value of the variational 
parameters a, ( are shown for both trial functions as a function of the interaction strength. 
Both the variational trial function of the LG condensate with vortex and without vortex are 
broadened under the influence of the repulsive interactions. 



III. VORTEX METASTABILITY BARRIERS 

A Laguerre-Gaussian condensate with a vortex can be metastable if there exists an energy 
barrier separating this state from a state without a vortex. In this section we estimate the 



height of the energy barrier both by the method proposed by Benakli et al. [£(J for 2D traps, 
and by the method proposed by Fetter and co-workers M. 



A. Hydrodynamic and microscopic instability 

The hydrodynamical and the microscopic instabilities studied by Fetter and co-workers 
@ involve the displacement of the vortex core relative to the center of the trap. As discussed 
in the introduction, for a parabolically trapped condensate it is energetically favorable for the 
distance between the vortex core and the center of the trap to increase. Thus, if dissipation is 
present, the vortex condensate in the parabolic trap can decay into a non- vortex condensate 
through a migration of the vortex to the edge of the cloud. This effect is difficult to observe 
experimentally, since the vortex core can tilt and the image contrast between the vortex 



core and the cloud of atoms is reduced 27 



To estimate the height of the energy barrier separating the vortex condensate from the 
condensate without a vortex in our present case, we introduce a third variational function 
to represent the order parameter of a vortex LG-condensate where the vortex core is at a 
given distance R from the cylindrical symmetry axis of the Laguerre- Gauss trapping beam: 
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MR; P, 0, z) = Af 2 r(p, 9) exp{-a R (p - p R ) 2 /2 - ( R z 2 /2}e^ 9 \ 



(6) 



The function r(p, 9) gives the distance between the point {p, 9, 0} and the vortex core: 



and the function 0(p, 6) gives the angle between the line connecting the vortex core with 
the point {p, 9, 0} and the a>axis: 



The parameters {&r, Pr, (r} are determined variationally for every studied distance R be- 
tween the vortex core and the center of the trapping geometry (the axis of propagation of 
the LG beam), and A/2 is a normalization constant. The limit R — > retrieves the results for 
the vortex condensate studied earlier. In the limit R — > 00, the result tends to the result for 
the condensate without a vortex. Using the variational approach, we calculate the energy 
E(R) of the displaced vortex state for any intermediate R. If the energy E(R) decreases 
monotonously with increasing R, the vortex is unstable - in the presence of dissipation the 
vortex condensate will decay. If there exists a maximum energy E(R max ) for an R max differ- 
ent from zero, there is an energy barrier Carrier = E(R mSuX ) — E(R = 0) which will hinder 
the vortex decay. 

In Figure 4, the energy barrier -Ebarrier is shown as a function of iVa scat /aHO; the dimen- 
sionless measure of the interaction strength used in the previous section. The calculations 
were performed for a configuration so that A = 5 HQ, Wq = 20 a 2 H0 . We found that a 
metastability barrier exists (-^barrier > 0) for Na scat / clho > 1-6 ± 0.1. This means that vor- 
tices are metastable only if the effective interaction strength (controlled by the scattering 
length, the number of atoms and the confinement strength) is large enough: the metasta- 
bility is induced by interactions. Upon further increasing the interaction strength above the 
threshold value, the metastability barrier increases. The variational energy E(R) is shown 
as a function of R in the insets: once for an interaction strength such that the vortex is not 
stable, and once for a situation where the vortex is metastable. A reasonable estimate of 




(7) 




(8) 
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the experimental parameters required to realize this trapping geometry gives ano — 4 /im 
and Q = 50 Hz. This implies that for 87 Rb, the critical number of trapped atoms neces- 
sary to make vortices metastable is of the order of 10 3 . Since one can reasonably expect 
to trap more than 10 5 atoms in the trap, vortices created in a Bose gas in the confinement 
potential mentioned above, will be well into the metastable regime (Na scat / 'duo ~ 625 for 
N = 500000). 

B. Uniform transition to a non- vortex state 

The other method which we used to study the metastability of the vortices in Laguerre- 
Gauss condensates, was proposed by Benakli et al. |26j for 2D traps with a an axial hole, 
punctured by an off-resonance laser beam. In this method, a trial solution for the Gross- 
Pitaevskii equation is constructed as a superposition between the vortex state and the state 
without a vortex : 9, z) = C ip(p, 9, z) + Ciip v (p, 9, z) where | C 1 2 H- \C\\ 2 = 1. Adapting 
the coefficients Cq and C\ of this superposition, the function changes from that of a LG 
condensate with a vortex to that without a vortex. If the energy of the intermediate states 
is found to be higher than the energy of the vortex state, this constitutes an energy barrier 
against uniform transition from vortex state to the state without a vortex, making the former 
metastable. 

Using the variationally optimized trial functions for the order parameter of a LG con- 
densate with a vortex and without a vortex, we found with the method of Benakli et al. 
that such a metastability barrier exists for Na sca t/ o-ho > 1.9 ± 0.1. The dependence of the 
energy barrier on Na scat / 'a ho is shown in Figure 5. In the inset of Figure 5, the variational 
energy of the trial function \1/ is shown as a function of | C 1 2 for an interaction strength 
below (a) and above (b) the threshold value for metastable vortices. The calculations were 
performed for a configuration so that Aq = 5 Ml, W$ = 20 clho- Note that the method of 
ref. [^6| followed here does not describe the decay of vortices mediated by the excitations, for 
example by nucleation of vortex loops. However, for a two-dimensional system punctured 
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by an off resonance laser beam, Benakli et al. [26 could show that excitation- mediated de 



cay of vortices only becomes appreciable for vortices in higher angular momentum states. 



Furthermore, for a purely one-dimensional toroidal system, Kagan et al. |28| showed that 
for a weakly interacting system at low temperature, the relaxation time for the decay the 
superfluid persistent current due to phonons is strongly suppressed. 

Both the method of Fetter and the method of Benakli et al. |26| give qualitatively the 



same result, namely that above a critical strength of interaction, expressed by Na sca t/aHO, 
the vortex state of the LG condensate is stabilized by a metastability barrier. Quantitative 
differences remain between the two methods, which is to be expected since they describe 
different possible mechanisms of vortex decay. The general conclusion - the stabilization of 
the vortex state due to the toroidal geometry - is further supported by a calculation for a 



vortex condensate in a mexican hat potential [29]. The metastability barrier calculated as 



by the method of ref. JF| is smaller than the barrier calculated used the method of ref. [26 



for Na scat /a HO > 3.5, which indicates that the latter method 0] overestimates the real 
metastability barrier. 

IV. INTERFERENCE AND DETECTION OF VORTICITY 

Phase coherence, a property necessary for the existence of vortices, was demonstrated 
experimentally using interference experiments f30fl . Two parabolically trapped condensates, 
displaced by a given distance, are allowed to expand freely. In the region where the two 
expanding condensates overlap, an interference pattern is observed. This type of experiment 



has also been proposed to observe vorticity in parabolically trapped condensates |17|,[18| : the 
presence of a vortex leads to an observable edge dislocation in the pattern of otherwise parallel 
interference fringes. 

Consider a parabolically trapped condensate positioned on the symmetry axis of the 
LG condensate. The LG condensate forms a cylindrical shell surrounding the paraboli- 
cally trapped condensate in the center. As the trapping potentials are switched off, both 
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condensates will expand and exhibit an interference pattern in the region of overlap. This 
interference pattern is the subject of the present section. 

To find the function representing the order parameter at a given time t after the start 
of the free expansion, the original function tp v (p, 9,z;t = 0) is expanded in free particle 
eigenfunctions. These eigenfunctions acquire a phase factor as time elapses, so that ip v 
at a time t is given by the resummed expansion with the 'time evolved' eigenfunctions. 
More explicitly, for the LG condensate with a vortex, first the coefficients of the plane wave 
expansion are evaluated: 

^ v (p,9,z;t = 0) = j c (k)j£^dk, (9) 

fu a u \ fir a ,exp{-ik p pcos(6 - <p) - ik z z} 
c{k p , (f), k g ) = J ip v {p, 6, z) (2tt) 3 / 2 ^ ' 



exp{-fc 2 2 /2C} 



x e l ^ l2) J dp p 2 .h{k p p) exp{-a(p - p ) 2 /2}. (11) 



o 



In this expression J±(x) is the Bessel function of first order of the first kind and the wave 
number k is expressed in cylindrical coordinates {k p ,<f),k z }. As time elapses, the free par- 
ticles eigenfunctions acquire a phase factor exp{ihk 2 t/(2m)}. The function at time t after 
the start of the free expansion is found by resumming the eigenfunctions at time t : 

/ e ik.r+ihk 2 t/(2m) 
c(k) ^j- 2 dk (12) 

oxp < - —. — - \ x e ie 



1 + <t [ 2(1 + i(t) 

oo 

x jdk p (p0 2 ^^e^'--)V 2exp {,^±|^} . (13) 
o 

The time evolution for the free expansion of a condensate prepared in a parabolic confinement 
is derived analogously |T7[ . The total measured density generated by the two condensates 



is then given by \ip v (p,9, z;t) + ip(p, 9, z;t)\ 2 - remember that each condensate function is 
normalized to the number of particles in the given condensate. Several time frames of the 
resulting evolution of the density of the expanding condensates are shown in Figure 3. The 
frames in Figure 3 show a cross-section of the density along the xy-plane. 
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Immediately after the traps are switched off, the density is that of a cylindrical, Laguerre- 
Gaussian condensate with a parabolically trapped condensate in the middle. As time goes 
by, both condensates expand: the parabolically trapped condensate expands radially, and 
the cylindrical shell of the LG condensate broadens. As the expanding condensates start 
to overlap the fringe pattern appears. If the LG condensate does not contain a quantum 
of superffuid circulation (i.e. no vortex), the interference pattern consists of a series of 
concentric circles with linearly increasing radius. If however the LG condensate does contain 
a vortex, the interference pattern is an Archimedean spiral. At higher vorticity, the number 
of arms in the Archimedean spiral equals the number of vortex quanta in the LG condensate. 

Figure 3 shows that, as a function of time, the spiral interference pattern rotates around 
the cylindrical symmetry axis of the trapped condensates, with a frequency of the order 
of the frequencies characterizing the parabolic approximation to the trapping potentials. 
Furthermore, as can be seen from Figure 3, the distance between two successive windings of 
the spiral increases as time increases. Hence, to detect the spiral interference pattern, the 
density has to be measured on a time scale shorter than that given by the inverse of the 
characteristic trapping frequencies. 



V. CONCLUSIONS 

In conclusion, we have studied the properties of vortices in a condensate in an optical 
trap generated by a laser beam in a Laguerre-Gaussian mode and proposed a method to 
detect these vortices. The energy of the vortex state was calculated variationally in a mean- 
field framework, and a threshold interaction strength was found beyond which there exists 
a metastability barrier stabilizing the LG vortex state against a transition to a non-vortex 
state. Both the hydrodynamic instability |7|] and the uniform transition to a non-vortex state 
26| were considered in the investigation of the vortex stability. Below the critical interaction 



strength the metastability barrier vanishes, irrespective of the fact that the present confined 
geometry allows for the vortex core to be positioned in a region where the condensate is 
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nonzero. 

The time evolution under free expansion of a LG condensate was derived, and sub- 
sequently used to calculate the interference pattern which arises from the overlap of an 
expanding LG condensate with an expanding parabolically trapped condensate in its center. 
When superfluid circulation is present, the interference pattern consists of an Archimedean 
spiral, clearly distinct from the series of concentric cylindrical fringes which arise in the 
interference pattern when no vorticity is present. The observation of such a spiral interfer- 
ence pattern would constitute clear evidence for vorticity in a spin-polarized Bose-Einstein 
condensates. 
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FIGURE CAPTIONS 



Figure 1 : The Laguerre-Gauss (LG) trap is an optical dipole trap consisting of a red- 
detuned laser beam. A typical intensity profile of the laser beam in the LG-mode {0, 1} is 
shown in this figure, in a cross section through the axis of propagation of the beam (the 
z-axis). The atoms will feel the optical dipole force attracting them to the (toroidal) region 
of highest intensity of the laser beam.v^O = 4. 472 14 



Figure 2 A typical surface of constant intensity of the laser beam in the Laguerre-Gauss 
(LG) mode {0, 1} is shown. The atoms in the red-detuned LG beam will be attracted to 
a toroidal or cylindrical region such as that within the surface of Figure 2. An additional 



parabolic magnetic confinement potential can be added in the ^-direction [ |l5| . |lq| 



Figure 3 : The variational result for the energy of the condensate in the Laguerre- 
Gaussian optical trap is shown as a function of the interaction strength. The full curve 
shows the energy of the condensate without a vortex, the dashed curve shows the energy 
with a vortex. The inset depicts the results for the variation parameters 00,01 in the trial 
function for the condensate with vortex (dashed line) and without (full line), as a function of 
interaction strength. The trapping (beam) parameters were chosen as follows: Aq = 5 HQ, 
Wq = 20 a 2 HQ . In all graphs of this Figure, energies are expressed in units Ml and lengths 
in units dfjo = 

Figure 4 : The energy barrier for removing a vortex from a condensate in a Laguerre- 
Gaussian optical trap, is depicted as a function of interaction strength. This metastability 
barrier ^barrier was calculated using the formalism of Fetter and co-workers [0. In this 
formalism, -Ebamer is found by deriving the energy E(R) as a function of the distance R 
between the vortex core and the center of the trapping potential, here the axis of propagation 
of the laser beam. This is illustrated in the insets. For interaction strengths lower than a 
critical value Na scat / duo < 1-6 ±0.1 the barrier vanishes and vortices in Laguerre-Gaussian 
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condensates are unstable with respect to the non-vortex state. For interaction strengths 
above this critical value, a metastability barrier exists. 



Figure 5 : The energy barrier for a uniform transition from a vortex condensate to a 
condensate without a vortex, in a Laguerre- Gaussian optical trap, is depicted as a function 
of interaction strength. The energy barrier for this mechanism was calculated by the method 



of Benakli et al. [Pq| . In the inset, the energy per particle of the variational function \I> = 
Coip + C\ip v is given as a function of | Co | 2 , for two different interaction strengths: one below 
(a) and one above (b) the threshold for metastability of vortices, and the relation to the 
energy barrier for a uniform transition is shown. 

Figure 6 : Several time frames in the evolution of freely expanding, overlapping con- 
densates are shown. The gray scale represents the density in a cross-section orthogonal to 
the symmetry axis of the trap (the direction of propagation of the Laguerre-Gaussian laser 
beam), with black being the maximum density and white the minimum density. The initial 
situation is depicted in the top left panel: a parabolically trapped condensate in the center 
is surrounded by a cylindrical shaped Laguerre-Gaussian condensate containing a vortex. 
The time evolution of the density in the region within the rectangle in the top left panel is 
shown in more detail in the subsequent panels. When the confinement of both condensates 
is switched off, they expand and overlap - the density in the gap between both condensates 
increases - and a spiral interference pattern is formed. Units are chosen so that lengths are 
expressed in a HO = \Jh/{mVl) and the time is in units of 1 / O. 
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